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We introduce a convenient parametrization of dark energy models that is general enough to 
include several modified gravity models and generalized forms of dark energy. In particular we 
take into account the linear perturbation growth factor, the anisotropic stress and the modified 
Poisson equation. We discuss the sensitivity of large scale weak lensing surveys like the proposed 
DUNE satellite to these parameters. We find that a large-scale weak-lensing tomographic survey 
is able to easily distinguish the Dvali-Gabadadze-Porrati model from ACDM and to determine the 
O^l ■ perturbation growth index to an absolute error of 0.02 — 0.03. 
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INTRODUCTION 



The observed late-time accelerated expansion (e.g. [lj) of the Universe has opened a Pandora's box in cosmology. 
Although a great number of models have emerged from the box, none of them has provided a satisfactory explanation 
of the observations ([2] and references therein). 

In a landscape of hardly compelling theories, probably the most efficient way to proceed is to exploit present and 
future data in search of signatures of unexpected phenomena that may signal new physical effects. In this way it is 
possible that we will be able to distinguish, say, a cosmological constant from dynamical dark energy (DE) or the 
latter from some form of modified gravity. 

To this end, an important task is to provide observational groups with simple measurable parameters that may 
be linked to interesting physics. In this paper we will investigate the extent to which additional parameters can be 
used to detect signatures of new cosmology in future surveys, with particular emphasis on weak lensing (WL). These 
additional parameters are connected to the growth of linear perturbations, to the anisotropic stress (defined as the 
difference bewteen the Newtonian potentials in longitudinal gauge) and to deviations from the Poisson equation for 
' matter. All these parameters find a simple motivation in current models of modified gravity, from extradimensional 
ones [ToJ to scalar-tensor theories; it is clear however that their introduction is not limited to these cases and can 
£SJ ' in fact account also for other phenomena, for instance clustering in the DE component. We also discuss which 
of the parameters arise naturally in which context. We then evaluate the sensitivity of WL experiments for two 
cases: a phenomenological one, in which the parameterization is chosen mainly on the grounds of simplicity and in 
I analogy to some specific models; and a more physically motivated case, namely the Dvali-Gabadadze-Porrati (DGP) 
<^ • [6J] extra-dimensional model. 

We focus on WL for two reasons: first, contrary to eg supernovae or baryon-oscillation tests, WL makes use of both 
background and linear perturbation dynamics, allowing to break degeneracies that arise at the background level (this 
is particularly important for testing modified gravity); second, several groups are planning or proposing large WL 
experiments in the next decade (e.g. 0,0, 0]) that will reach the sensitivity to test cosmology at unprecedented depth 
and it is therefore important to optimize the science return of these proposals. We will therefore produce Fisher matrix 
confidence regions for the relevant parameters for surveys like those proposed by the DUNE @] or JDEM/SNAP 0,13] 
collaborations. 

The effect of the anisotropic stress has been discussed several times in DE literature. Refs. [131 and [2l| discuss 
it in the context of scalar-tensor theories and evaluated the effect on the shear spectrum; [36] discusses its role for a 
coupled Gauss-Bonnet theory; [13] showed that it is essential to mimic the DGP growth-rate with a fluid dark energy 
model; finally, during the final stages of this project another paper discussing the anisotropic stress in dark energy 
models appeared Similarly, the parametrization of the growth factor in DE models has been often investigated in 
the past PJ, [Bi- 



ll. DEFINING THE DARK SIDE 



In this section we discuss our parametrisation of the dark sector. We separate it into two parts, firstly the parametri- 
sation of the background quantities relevant for a perfectly homogeneous and isotropic universe, and secondly the 
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additional parameters which describe small deviations from this idealised state, at the level of first order perturbation 
theory. 

We assume throughout that the spatial curvature vanishes, mostly for simplicity. We note however that especially 
in the context of general dark energy models curvature is not yet well constrained [§] . 



A. Parametrisation of the expansion history 

In order to characterise the evolution of the Universe at the background level, we need to provide a parametrisation 
of the expansion history, H{z). To this end we write it as 

H{z) 2 = H 2 [Q m (l + zf + (1 - Q m )f(z)] . (1) 

We define O m by the requirement that during matter domination, i.e. at redshifts of a few hundred, the Universe 
expands according to 

H(z~ 0(100)) = H a ^TL^{l + z)i . (2) 

We will call this part the (dark) matter, and the remaining part the dark energy. However, these are only effective 
quantities 0. It could be that we are dealing with a tracking dark energy that scales like matter at high redshifts 
[37L l38|. In this case, the scaling part of its energy density would be counted as matter. It is also possible that the 
observed acceleration of the expansion is due to a modification of gravity. In this case there is no dark energy at all 
and the (1 — fl m )f(z) part takes into account the gravity terms. In all cases, the fundamental quantity is H(z). 

The function f(z) describes the evolution of the energy density in the (effective) dark energy. The most widely 
used parametrization makes use of a Taylor series in a [39], 

w(a) = w + (1 - a)w a . (3) 

where a = 1/(1 + z) is the scale factor. For the above form of w there is an analytic expression of /: 

f(z;w ,w a ) = (l + z)3( 1 +-o+-») e xp|-3u )aT ^-| (4) 

However this is not necessarily a good fit to physical models or to the data j4jj; in the DGP case we will make use of 
the full expression for w(a). 



B. Parametrisation of the first order quantities 



We will limit ourselves to scalar perturbations. In Newtonian gauge we can then write the line element defining the 
metric as 

ds 2 = a(r) 2 [-(1 + 2i>)dr 2 + (1 - 2(j))dx 2 ] (5) 

so that it is characterised by two perturbative quantities, the scalar potentials ip and 4>. 

The perturbations in the dark fluids are characterised for example by their comoving density perturbations A = 
S + 3HaV/k 2 and their velocity perturbations V. Their evolution is sourced by the potentials 4> and ip and depends 
also on the pressure perturbations dp and the anisotropic stresses a of the fluids. For the (dark) matter we set Sp = 
(as in addition w m = this is a gauge-invariant statement) and a m — 0. In App. [A] we define our notation and review 
the perturbation formalism in more detail. 

In General Relativity the <f) potential is given by the algebraic relation 

-k 2 dp = AirGa 2 Y J P l ^ (6) 

i 

which is a generalisation of the Poisson equation of Newtonian gravity (the factor —k 2 is the spatial Laplacian). We 
see that in general all fluids with non-zero perturbations will contribute to it. Since we have characterised the split 
into dark matter and dark energy at the background level, we cannot demand in addition that Ade "C A m . At 
the fluid level, the evolution of A DE is influenced by a combination of the pressure perturbation and the anisotropic 
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stress of the dark energy. However, the pressure perturbation of the dark energy is only very indirectly related to 
observables through the Einstein equations. For this reason we rewrite Eq. ^ as 

k 2 4> = -4irGa 2 Qp m A m . (7) 

where G is the gravitational constant measured today in the solar system. Here Q(k, a) is a phenomenological quantity 
that, in general relativity (GR), is due to the contributions of the non-matter fluids (and in this case depends on their 
Sp and a). But it is more general, as it can describe a change of the gravitational constant G due to a modification of 
gravity (see DGP example below). It could even be apparent: If there is non-clustering early quintessence contributing 
to the expansion rate after last scattering then we added its contribution to the total energy density during that period 
wrongly to the dark matter, through the definition of Q m , In this case we will observe less clustering than expected, 
and we need to be able to model this aspect. This is the role of Q{k, a). 

For the dark energy we need to admit an arbitrary anisotropic stress a, and we use it to parametrise ip as 

^ = [1 + I7(M)]0. (8) 

At present there is no sign for a non- vanishing anisotropic stress beyond that generated by the free-streaming of photons 
and neutrinos. However, it is expected to be non-zero in the case of topological defects [ll[ or very generically for 
modified gravity models [T3|. In most cases, -q(a — > 0) — > in order to recover the behaviour of the standard model, 
however in some models like scalar-tensor theories this may not be the case. 

If Q ^ 1 or r] ^ 0, then we need to take into account the modified growth of linear perturbations. Defining the 
logarithmic derivative 

dlogA m 

— = m[a) (9) 

ci log a 

it has been shown several times [16 , [13, H, 0, Ell that for several DE models (including dynamical dark energy and 
DGP) a good approximation can be obtained assuming 

m{a) = tt m {ay (10) 

where 7 is a constant that depends on the specific model. 

Although for an analysis of actual data it may be preferable to use the parameters {Q,ri}, we concentrate in this 
paper on the use of weak lensing to distinguish between different models. As we will see in the following section, the 
quantities that enter the weak lensing calculation are the growth index 7 as well as the parameter combination 

£ = Q(l + ?/ /2). (11) 

Weak lensing will therefore most directly constrain these parameters, so that we will use the set {7, £} for the 
constraints. Of course we should really think of both as functions of Q and 77. An additional benefit of using 7 is that 
as anticipated it is relatively easy to parametrise since in most models it is reasonable to just take 7 to be a constant 
(see section [IIIB|) , For £ the situation is however not so clear-cut. First, a simple constant £ would be completely 
degenerate with the overall amplitude of the linear matter power spectrum, so would be effectively unobservable per 
se with WL experiments. Then, even simple models like scalar-tensor theories predict a complicate time-dependence 
of E so that it is not obvious which parametrization is more useful (while for DGP it turns out that E = 1 just as in 
GR, see below). So lacking a better motivation we start with a very simple possibility, namely that E starts at early 
times as in GR (i.e. E = 1) and then deviates progressively more as time goes by: 

E(a) = l + E a. (12) 

This choice is more general than it seems. If we assume that E(a) is linear and equal to unity at some arbitrary a\, 
then we have E(a) = 1 + Ei(ai — a) = (1 + Eiai)[l — aEi/(l + Eiai)]. But since an overall factor can be absorbed 
into the spectrum normalization, one has in fact E(a) = 1 — E^a where £'1 = Ei/(1 + Eiai). Therefore, from the 
error on So in lfT2|) one can derive easily the error on the slope Ei at any point a±. 

As a second possibility we also investigate a piece-wise constant function E(a) with three different values £1,2,3 in 
three redshift bins. Here again, due to the degeneracy with as, we fix £1 = 1. 

In this way we have parametrised the three geometric quantities, H, <f> and ip, with the following parameters: Q rn , £0 
(or £2,3) and 7, plus those that enter the effective dark energy equation of state. These represent all scalar degrees of 
freedom to first order in cosmological perturbation theory. The parameterisations themselves are clearly not general, 
as we have replaced three functions by six or seven numbers, and as we have not provided for a dependence on k at all, 
but they work well for ACDM, Quintessence models, DGP and to a more limited extent for scalar-tensor models. We 
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emphasize that this set is optimised for weak lensing forecasts. For the analysis of multiple experiments, one should 
use a more general parametrisation of Q and r\. If the dark energy can be represented as a fluid, then {Sp, a} would 
be another natural choice for the extra degrees of freedom, with rj = rj(a) and Q = Q(8p, a). While phenomenological 
quantities like Q, rj, £ and 7 are very useful for measuring the behaviour of the dark energy, we see them as a first 
step to uncovering the physical degrees of freedom. Once determined, they can guide us towards classes of theories - 
for example 77 = would rule out many models where GR is modified. 



III. OBSERVABLES 



A. Constraining the expansion history 

In order to constrain the expansion history we can measure either directly H (z) or else one of the distance measures. 
The main tools are 

• Luminosity distance: Probed by type- la supernova explosions (SN-Ia), this provides currently the main 
constraints on H(z) at low redshifts, z < 1.5. 

• Angular diameter distance: Measured through the tangential component of the Baryon Acoustic Oscillations 
(BAO) either via their imprint in the galaxy distribution [HS] or the cosmic microwave background radiation 
(CMB) [321 ] . One great advantage of this method is the low level of systematic uncertainties. The CMB provides 
one data point at z « 1100, while the galaxy BAO probe mostly the range z < 1.5. Using Lyman-a observations 
this could be extended to z « 3. 

• Direct probes of H(z): This can either be done through the radial component of the BAO in galaxies, or 
through the dipole of the luminosity distance [28]. 



B. Growth of matter perturbations 

In the standard ACDM model of cosmology, the dark matter perturbations on sub-horizon scales grow linearly with 
the scale factor a during matter domination. During radiation domination they grow logarithmically, and also at late 
times, when the dark energy starts to dominate, their growth is slowed. The growth factor g = A m /a is therefore 
expected to be constant at early times (but after matter-radiation equality) and to decrease at late times. In addition 
to this effect which is due to the expansion rate of the universe, there is also the possibility that fluctuations in the 
dark energy can change the gravitational potentials and so affect the dark matter clustering. 

In ACDM g can be approximated very well through 

g(a) = ex p{^ d\na(fl m (a)' 1 - 1)| (13) 

where 7 ~ 0.545. There are two ways that the growth rate can be changed with respect to ACDM: Firstly, a general 
w of the dark energy will lead to a different expansion rate, and so to a different Hubble drag. Secondly, if in the 
Poisson equation ([7j we have Q ^ 1 then this will also affect the growth rate of the dark matter, as will 77 7^ 0. We 
therefore expect that 7 is a function of w, 77 and Q. 

For standard quintessence models we have 77 = and Q « 1 on small scales as the scalar field does not cluster due 
to a sound speed c 2 s — 1. In this case 7 is only a weak function of w fis| : 

7(10) = 0.55 + 0.05(1 + w(z = 1)). (14) 

A constant 7 turns out to be an excellent approximation also for coupled dark energy models [35| and for modified 
gravity models A similar change in 7 can also be obtained in models where the effective dark energy clusters 

If we assume that ip is the dominant source of the dark matter clustering, then we see that its value is modified 
by (1 + rf)Q (while the 4> term sourcing perturbations in 6 m does not contain the contribution from the anisotropic 
stress). Following the discussion in [2(3] we set 

1 - ilm{a) 
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and derive the asymptotic growth index as 

= 3(1 -woo -A(Q,r))) ^ 
5 - 6woo 

In section [TV PI on the DGP model we show how these relations can be inverted to determine <f> and ip . 

Although these relations are useful to gain physical insight and for a rough idea of what to expect, for an actual 
data analysis one would specify Q and rj and then integrate the perturbation equations. 



C. Weak lensing 

Usually it is taken for granted that matter concentrations deflect light. However, the light does not feel the 
presence of matter, it feels the gravitational field generated by matter. In this paper we consider a scenario where the 
gravitational field has been modified, so that we have to be somewhat careful when deriving the lensing equations. 
Following [2l|| we exploit the fact that it is the lensing potential <£> = <p + 4> which describes the deviation of light rays 
in our scenario. In ACDM we can use the Poisson equation to replace the lensing potential with the dark matter 
perturbations (since the cosmological constant has no perturbations), 

fc2$ = 2 ^k Aro . (17) 
la 

However, in general this is not the complete contribution to the lensing potential. For our parametrisation we find 
instead that 

kH = Q(2 + r,)^I^A m = 2£^f^A m . (18) 
2a la 

Correspondingly, the total lensing effect is obtained by multiplying the usual equations by a factor £ = (1 + j]/2)Q. 
We also notice that a modification of the growth rate appears twice, once in the different behaviour of A m , and a 
second time in the non-trivial factor S. It is clear that we need to take both into account, using only one of the two 
would be inconsistent. 

An additional complication arises because some of the scales relevant for weak lensing are in the (mildly) non-linear 
regime of clustering. It is therefore necessary to map the linear power spectrum of <£> into a non-linear one. This is 
difficult even for pure dark matter, but in this case fitting formulas exist, and numerical tests with iV-body programs 
have been performed. For non-standard dark energy models or modified gravity the mapping is not known, and it 
probably depends on the details of the model. Although this may be used in the long run as an additional test, for 
now we are left with the question as to how to perform this mapping. In this paper we have decided to assume that 
SA m represents the effective clustering amplitude, and we use this expression as input to compute the non-linear 
power spectrum. 

One side effect of this prescription is that as anticipated any constant pre-factor of A m becomes degenerate with 
as (which is also a constant) and cannot be measured with weak lensing. In practice this means that if we use {7, £} 
as the fundamental parameters, then a constant contribution to £ cannot be measured. On the other hand if we use 
{Q, 77} as fundamental parameters then even a constant 77 affects 7 and the degeneracy with as is broken through this 
link. 



D. Other probes 



We have already mentioned one CMB observable, namely the peak position, which is sensitive mostly to the 
expansion rate of the universe and provides effectively a standard ruler. On large scales the CMB angular power 
spectrum is dominated by the integrated Sachs- Wolfe (ISW) effect, which is proportional to <fi + ip = $. The ISW 
effect therefore probes a similar quantity as weak lensing, but is sensitive to its time evolution. Additionally, it is 
mostly relevant at low £, which means that it is strongly affected by cosmic variance. This limits its statistical power 
in a fundamental way. 

From the perturbation equation for the matter velocity perturbation, Eq. (|A15|) we see that V is only sensitive to 
tp. Additionally, the peculiar velocities of galaxies are supposed to be a very good tracer of the dark matter velocity 
field [IB]- The peculiar velocities are therefore a direct probe of 1/) alone. This makes it in principle an excellent 
complement of weak lensing (which measures <f> + ip) and of the growth rate (measuring mostly <f> through the Poisson 
equation) but of course measuring reliably peculiar velocities to cosmological distances is still prohibitive. 
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The perturbations in the metric also affect distance measurements like e.g. the luminosity distance [27l. |29|. The 
fluctuations in the luminosity distance on small angular scales are a measure of both the peculiar motion of the 
supernovae as well as the lensing by intervening matter perturbations. As very large supernova data sets are expected 
in the future, this may turn into a promising additional probe of cj> and -0. 

IV. DARK ENERGY MODELS 

Let us now review some of the models among those presented in the literature that are amenable to our parametriza- 
tion. 



A. Lambda-CDM 

In ACDM the dark energy contribution to the energy momentum tensor is of the form T" = A.6". Prom this 
we see immediately that p — — p so that we have a constant w = — 1, Additionally there are no perturbations in 
a cosmological constant, which is compatible with the perturbation equations (| A8[) and (| A9[) as they decouple from 
the gravitational potentials for w = — 1. For this reason Q = 1. The absence of off-diagonal terms in the energy 
momentum tensor also shows that 77 = so that S = 1. The growth-rate of the matter perturbations is only affected 
by the accelerated expansion of the universe, with 7 w 6/11. 

B. Quintessence 

In quintessence models the dark energy is represented by a scalar field, so that w can now vary as a function of 
time, subject to the condition w > —1 if pq > 0. At the level of first order perturbations, quintessence is exactly 
equivalent to a fluid with (? s = 1 and no anisotropic stress (?/ = 0). For such a fluid one often relaxes the condition 
on w, so that w < — 1 becomes allowed (see e.g. [HI and papers cited therein for more details). The high sound 
speed suppresses clustering on sub-horizon scales so that \Q{k 3> Hq) — 1| <C 1. It is on these scales that weak lensing 
and galaxy surveys measure the matter power spectrum, so that E w 1, and the growth-rate is only slightly changed 
through the different expansion rate. From (fl4"ll we see that typically one has 7 = 0.54 -7- 0.57 for the range of w still 
allowed by the data. 

On very large scales clustering can be non-negligible, which affects for example the low I part of the CMB spectrum 
through the ISW effect (25|. For a DGP-like equation of state, we find Q(a = 1, k — Hq) < 1.1 by integrating the 
perturbation equations numerically. In general we expect that Q — 1 remains small below the sound horizon, while it 
can be non-negligible on larger scales. Thus, if we find that Q « 1 on small scales, but Q > 1 for scales larger than 
some critical scale, then this may be a hint for the existence of a sound horizon of the dark energy. Additionally, 
one important class of quintessence models exhibits a period of "tracking" at high redshift, during which the energy 
density in the dark energy stays constant relative to the one in the dark matter [37], [H| . If the dark energy makes up 
an important fraction of the total energy density during that period, O m as defined through Eq. ^ will be too large. 
In that case Q < 1 as the quintessence clusters less than the dark matter on small scales. 

C. A generic dark energy model 

A generic non-standard model with either a dark energy very different from quintessence or in which gravity is non- 
Einsteinian will presumably introduce modifications in the Poisson equation, the equation for dark matter growth 
and/or the anisotropic stress. Since these quantities are not independent, a generic modified gravity model should 
allow for at least an additional parameter for E and one for 7. As we anticipated we assume 7 constant and then 
either E = 1 + Eoa or a piece-wise constant E in three redshift bins. Let us denote these two "generic dark energy" 
models as GDE1 and GDE2. 

Since these parametrisations contain ACDM in their parameter space, which is the phenomenologically most suc- 
cessful model today, they are useful to characterise the sensitivity of experiments to non-standard dark energy models. 
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D. DGP 



An alternative approach to the late-time accelerated expansion of the universe modifies the geometry side of 
the Einstein equations, rather than the energy content. A well-known model of this kind is the Dvali-Gabadadze- 
Porrati model Q. This model is based on five-dimensional gravity, with matter and an additional four-dimensional 
Einstein-Hilbert action on a brane. This then modifies the evolution of the Universe, with one solution asymptotically 
approaching a de Sitter universe. Assuming spatial flatness, the Hubble parameter for this solution is given by 



H 2 



We can solve this quadratic equation for H to find 



H 



8ttG 




(19) 



(20) 



Since the matter stays on the brane, its conservation equation is four-dimensional so that p m oc a -3 . Considering an 
effective dark energy component with p e g = 3H/ (87rGr c ) which leads to the DGP expansion history, one can use the 
conservation equation p c g + 3H(1 + w c s)p c s = to define an effective equation of state with |l2j 



w eS (a) 



Om - 1 - y/(l - On) 2 + 4» m /q 3 

2v , (l-On) 2 +4f! m / a 3 



(21) 



For Q m = 0.3 we find wq » —0.77 and w a « 0.3. For 77 we turn to [13, [la] which have computed the perturbations on 
small scales and found 



k z (t> = -AnGa 2 ( 1 - — I p m A 



where the parameter {3 is defined as: 

13 = 1- 2Hr c \ 1 



3ft 



iirGa 2 [ 1 + — ) p m A 



H 
3IP 



2(Hr c 
2Hr r - 



1 + 2i?r c w DE 



(22) 
(23) 

(24) 



With the parametrisation Eq. ((SJ we find that 



V 



Q m (a) 



3/3-1 2 + n m ( a y 



(25) 



which means that its value today is r/(a = 1) w —0.44 for Cl m = 0.3. We plot its evolution as a function of a in Fig. [TJ 
We see that it vanishes at high redshift, when the modifications of gravity are negligible. 
For the growth function we turn to [2(j]. They find 



7 + 5q Tt (a) + 7r; 2 » + 3» 3 „( a ) 
[1 + n^(o)][ll + 50 m (a)] 



(26) 



Fig. [2] compares this formula with the numerical result, and we see that it works very well. On average in the range 
z G (0, 3) we can use 7 w 0.68. 

Q is just the non-trivial pre-factor in Eq. (|22| . 



= 1 



1 

3^5' 



(27) 



it is plotted in Figure [3 We find that for DGP S = 1 so that the deviation of light rays due to a given mass is the 
same as in GR. The weak lensing results therefore depend only on 7 which modifies the growth of A m , but not on Q 
and 77 separately. In this argument we also used that the gravitational constant G C av measured by a Cavendish type 
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Figure 1: The r\ parameter of the anisotropic stress as a function of a for Q m = 0.3. The black solid line shows the actual value 
while the red dashed curve shows the recovered rj using the fitting formula for 7. 




Figure 2: The growth parameter 7 of DGP, comparison between the fitting formula (red dashed curve) and the numerical result 
(black solid line) for Q m = 0.3. 



experiment is just the bare constant G as the force-law modifications of DGP are screened for A m 3> 1 This is 
in contrast to the situation for scalar-tensor theories. 

It is maybe instructive to illustrate with a short "Gedanken" experiment how we could recover 4> and i/i from data 
if our universe was described by a DGP-like model. We assume that w(a) has been measured by e.g. supernovae 
and/or BAO. Suppose now that looking at the growth-history of the matter power spectrum (available potentially as 
a by-product from a BAO survey or from a dedicated galaxy survey) we notice that 7 is outside the range allowed for 
quintessence-like models. We therefore have to assume that the dark energy is non-standard. We now need to invert 
the computation of 7 to get the link with 4>, through the definition of Q(a), 



-4irGa 2 Q(a)p m A r 



(28) 
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Figure 3: The Q parameter of DGP, fi m = 0.3. The black solid curve shows the exact value while the red dashed line is the 
result recovered with the fitting formula for 7. 



Using the relations in section IIII Bl we find 



(1 + ri(a))Q(a) - 1 = (1 - n m (a)) ( 1 - w(a) - ^J^ 7 (a) 



Lensing on the other hand gives 

k 2 ((j> +ip) = k 2 4>{2 + V ) = -(2£)47rGa 2 p m A m . 
For DGP the lensing signal is precisely the one expected naively for A TO , so that 

E = l<^?7= |- - 2 



(29) 



(30) 



(31) 



We find that (l + r?)Q = 2 — Q, allowing us to recover Q and rj separately. They are shown as the dashed curves in 
Figs. [Hand [H Once we know 77, 7 and Q we can compute <f> and ip. 

We find that the accuracy of the fitting formula for 7(Q, rf) is quite good, and certainly sufficient for current 
experiments, and for distinguishing DGP from Quintessence at the perturbation level. We will see later that next- 
generation weak lensing experiment can reach a precision where the differences are important. At that point one may 
need to numerically integrate the perturbation equations to compute 7. 



E. ADGP 



We find it instructive to introduce a simple variant of DGP that we denote as ADGP, namely a DGP model 
which includes a cosmological constant. In this way we can interpolate between DGP proper (Oa = 0) and ACDM 
= 1 — Cl m , ie fl rc = 0). Notice that we are still taking the self-accelerating branch of DGP, different from e.g. 
[1J|. In flat space, the Hubble parameter is then given by 



8ttG 



(Pm + PA) 



with the self-accelerating solution now being 



H=^- + \^ + — ( Pm + PA) = H l 



(32) 



(33) 
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In flat space we additionally have that Q Te = \ (1 — fl m — ft A ) 2 . We also notice that we can define an overall effective 
dark energy fluid through 



Peff = PA 



3H 

8ttGV c 



(34) 



from which we can derive an effective equation of state. Concerning the perturbations, it was shown in [14[ that the 
DGP force laws lf22|) and ([23)1 are not changed through the addition of a brane cosmological constant if we write (3 as 



(3=1- 2r c H ^1 

which depends on the value of A. For our key quantities we find 

1 - Q m (a) + S1a (a) 



H 
3H 2 



w(a) 
(3(a) 
T)(a) 
Q(a) 



[1 - Q m (a)] [1 + fl m (a) + Q A (a)] 

[1 + n m (a) + Qa (a)} 2 - 2Q m (a) 
[1 - il m (a) - n A (a)] [1 + Q m (a) + n A (a)] 

1 - [Q m (a) + r> A (a)] 2 

30 m (a) - [1 + Q m (a) + fi A (a)] [2 + VL m (a) + (a)] 

2 1 + n A (a) - 2Q m (a) + (1 + gm (a) + n A (a)) 2 

3 (1 + fi m (a) + Q A (a)) 2 - 2ft m (a) 



(35) 

(36) 
(37) 
(38) 



E(a) = 1 

These reduce to the ones of DGP and ACDM in the respective limits. The growth factor can be approximated by 



(39) 
(40) 



7(a) 



(i + n m + n A f - 2n m 


[7 - 4f> m + 5r> A - 3«i _ 3n m n A ] + 2n m (i - n m - n A ) 




(i + n m + n A ) 2 — 2fi m 


[(ii + 50 ro ) (i - n m ) + (ii - 50 ro ) n A ] 



(41) 



where for ease of notation we suppressed the explicit dependence of the fls on a 



F. Scalar-tensor theories 

For completeness we also give a brief overview of the relevant quantities in scalar-tensor theories (see eg. [HI, 01). 
For a model characterized by the Lagrangian 

L = F(tp)R - if.^ip ; " - 2V{ip) + \^G*L matter (42) 

(where F(ip) is the coupling function, that we assume to be normalized to unity today and G* is the bare gravitational 
constant) the relation between the metric potentials is 

1> = <t> ~ y&<P (43) 

where F' = dF/dip. It turns out that in the linear sub-horizon limit the functions tjj,<p,6(f obey three Poisson-like 
equations: 

k 2 4> 
k 2 5ip 

Then from |28|) we derive 



G* 2 . 2(F + F' 2 ) 

= — 47T a P m lA m 7T 

F 2F + 3F' 2 

G* 2 . 2(F + 2F' 2 ) 
= -At:— a p m A m 2p + 3pl2 

G* 2 . 2FF' 

= 47r — a p m A m - 

F 2F + 3F' 2 



(44) 
(45) 
(46) 



Q 



G* 2(F + F' 2 ) 
FG cav n 2F + 3F'2 



(47) 
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Model 


growth index 


E(a) 


new param. 


fid. values 


GDE1 


7 = const 


E = 1 + E a 


7, So 


(0.55,0) 


GDE2 


7 = const 


E(i-th z-bin) = E* 


7, S 2 ,E 3 


(0.55,1,1) 


DGP 


7 = const 


E = 1 


7 


0.68 


ADGP 


7 (a) 


E = 1 




or 0.7 



Table I: The DE models considered in this paper. 



where G cav ^ is the presently measured value of the gravitational constant in a Cavendish-like experiment. If the 
equation lj45|) can be assumed to hold in the highly non-linear laboratory environment then one would define 

_G* 2(Fo + 2iff) 



Moreover, we obtain the anisotropic stress 



Finally, we derive 



ib -6 F 12 



Q ( 1 + l) = FG 



s = 

(notice that our result differs from [ll|). It is clear then that depending on F our simple phenomenological parametriza- 
tion may be acceptable or fail completely. Moreover we find that the usual growth fit lfl3|) is not a very good approx- 
imation since during the matter era the growth is faster than in a ACDM model. The analysis of specific examples of 
scalar-tensor models is left to future work. 



V. FORECASTS FOR WEAK LENSING LARGE-SCALE SURVEYS 



We finally are in position to derive the sensitivity of typical next-generation tomographic weak lensing surveys to 
the non-standard parameters introduced above, expanding over recent papers like Refs. and [ifj. In particular, 
we study a survey patterned according to the specifications in Ref. [45|], which dealt with the standard model. In 
Appendix B we give the full convergence power spectrum as a function of rj and Q and the relevant Fisher matrix 
equations. 

Let us then consider a survey characterized by the sky fraction / s ky, the mean redshift z mcan « 1.412zo and the 
number sources per arcmin 2 , d. When not otherwise specified we assume z mcan = 0.9 and d = 35 as our benchmark 
survey: these values are well within the range considered for the DUNE satellite proposal. The derived errors scale 

1 /2 

clearly as fj^ so it is easy to rescale our results to different sky fractions. We assume that the photo- z error obeys 
a normal distribution with variance a z = 0.05. We choose to bin the distribution out to z = 3 into five equal-galaxy- 
number bins (or three for the model with a piece- wise constant £(z)). For the linear matter power spectrum we adopt 
the fit by Eisenstein & Hu [52] (with no massive neutrinos and also neglecting any change of the shape of the spectrum 
for small deviations around w = —1). For the non-linear correction we use the halo model by Smith et al. [43]. We 
consider the range 10 < i < 20000 since we find that both smaller t and larger ts do not contribute significantly. 

We begin the discussion with the generic dark energy models GDE1 and GDE2. The parameter set (with the 
fiducial values inside square brackets) is therefore 

Pa = {u> m = n m h 2 [0.U7],u} b = fi 6 /i 2 [0.02205],n s [l],r!„ l [0.3],wo[-0.95],Wa[0],7[0.55], ( T 8 [0.8]} (51) 

while for S we assume as fiducial values either So = (GDE1) or Si, 2, 3 = 1 (GDE2). 

First we study how the estimate of wo,w p (projection of wo,w a on the pivot point a p defined as the epoch at 
which the errors decorrelate) is affected by fixing the other parameters. In Fig. [4] we show the FOM defined as 
l/[a(wo) ■ a(w p )] first when all the parameters are fixed to their fiducial value (first bar) and then successively 
marginalizing over the parameter indicated in the label and over all those on the left (eg the fourth column represents 
the marginalization over uj m ,ui b ,n s ). This shows that the WL method would benefit most from complementary 



12 



n m h 2 fl b h 2 n fl m cr 8 y Z 

4000 
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2000 
1000 





n m h 2 n b h 2 n H m <r 8 7 S 

Figure 4: FOM for wo,w p vs. marginalized parameters of the model GDE1. 

experiments that determine Q mi il m h 2 . On the other hand, there is not much loss in marginalizing over the two 
non-standard parameters So, 7. 

In Fig. [5] we show the confidence regions for So, 7 . Errors of the order or 0.1 for So and 0.3 for 7 are reachable 
already with the benchmark survey. In Fig. [6] we show the FOM (wo,w p ) varying the depth z m ean ~ 1.412zo and 
the density d of sources per arcmin 2 (full marginalization) . If we set as a convenient target a FOM equal to 1000 (for 
instance, an error of 0.01 for wq and 0.1 for w p ) then we see that our benchmark survey remains a little below the 
target (we obtain a(wo) = 0.018 and a(w p ) = 0.088), which would require at least d = 50 or a deeper survey. 

In Fig. [7] we show the FOM (7, S ) again varying z mean , d (full marginalization). Here we set as target a FOM of 
5000, obtained for instance with an error 0.02 on 7 and 0.01 on S ; it turns out that the target can be reached with 
the benchmark survey. 

For the model GDE2 we divide the survey into three equal-galaxy-numbers z-bins and choose a S(z^) = S.; piece- 
wise constant in the three bins. Fixing Si = 1, we are left with two free parameters S 2 ,3. In Fig. [8] we show the 
confidence regions; we see that WL surveys could set stringent limits on the deviation of S from the GR fiducial value. 

We can now focus our attention to the DGP model. As anticipated, in order to investigate the ability of WL studies 
to distinguish the DGP model from ACDM, we consider two cases. First, we assume a standard DGP model with 
w{z) given by Eq. (|2"Tjl . In this case the model also determines the function 7(2). For S! m w 0.3 one has an almost 
constant 7 in the range z 6 (0, 3) with an average value 7 « 0.68. Instead of using the full equation for j(z) we prefer 
to leave 7 as a free constant parameter in order to compare directly with a standard gravity DE model with the same 
w(z) and the standard value 7 « 0.55. In Fig. [9] we show the confidence regions around the DGP fiducial model; our 
benchmark surveys seems well capable of differentiating DGP from ACDM. 

Then, we consider the ADGP model, in which the "matter" content is in fact matter plus a cosmological constant, 
so that in the limit of fl\ = one recovers DGP, while when Q\ = 1 — O m one falls back into pure ACDM. In Fig.fTOl 
we see again that our benchmark survey will be able to distinguish between the two extreme cases with a very high 
confidence. In Fig. Q1] we display the weak lensing spectrum for ACDM in the 5th z-hin with the noise due to the 
intrinsic ellipticity and for comparison the DGP spectrum. We see that the DGP spectrum is well outside the noise 
at low ^'s. 



VI. CONCLUSIONS 

In this paper we propose a general parametrisation of both dark energy and modified gravity models up to the linear 
perturbations. Apart from parametrising the Hubble parameter H(z) with an effective equation of state, we use the 
growth index 7 and the effective modification of the lensing potential S. We discuss the relation of these quantities to 
the anisotropic stress (parametrised through 77 or a) and the modification of the Poisson equation (given in terms of a 
parameter Q). We then show how these parameters appear in different experimental setups, concentrating specifically 
on the case of weak lensing. We also give explicit expressions for the parameters for a range of models like ACDM, 
Quintessence, the DGP model and scalar-tensor theories. We identify a few signatures that could point to specific 
theories: The detection of a significant anisotropic stress would favour modified-gravity like theories, while strong 
upper limits on \r]\ could rule out many such models. A significant deviation from Q — 1 on large scales only could 
point to a finite sound speed of the dark energy. The hope is that we can eventually use such clues to understand the 
physical nature of the phenomenon underlying the accelerated expansion of the Universe. 
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Figure 5: Confidence regions at 68% for the benchmark survey z mea n — 0.9, d — 35 (outer contour) and for d = 50, 75 (inner 
contours). 
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Figure 6: FOM (wo,w p ) vs. z mean and d (galaxies per arcmin 2 ) in GDE1. The grey area represents a convenient target. 
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Figure 7: FOM (7, So) vs. z mean and d (galaxies per arcmin 2 ) in GDE1. The target is within the grey area. 
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0.940.960.98 1 1.021.041.06 

Figure 8: Confidence regions at 68% for the benchmark survey z mea n — 0.9, d = 35 (outer contour) and for d = 50, 75 (inner 
contours) in the model GDE2. 
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0.55 
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Figure 9: Confidence regions at 68% for the benchmark Survey Zmean 
contours) for DGP. The dotted line represents the ACDM value. 



0.9, d = 35 (outer contour) and for d = 50, 75 (inner 



We use our parametrisation to provide forecasts for weak lensing satellite experiments (having in mind a setup 
similar to DUNE) on how well they will be able to constrain dark energy and modified gravity models. We find that a 
DUNE-like survey will be able to constrain the growth index with an error that varies from 0.015 to 0.036 depending 
on the model (see Table Hlf. This is sufficient to rule out a model like DGP at more than 7 standard deviations based 
on the perturbations. Table ITTI also shows that the parameter £ can be strongly constrained. This demonstrates that 
weak lensing will evolve in the next decade into a very powerful probe of the dark energy phenomenon, with the 
potential to deliver insight into the physics behind the accelerated expansion of the universe through constraints on 
the dark energy perturbations and on modified gravity. 
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Figure 10: Confidence regions at 68% for the benchmark = 0.9, d = 35 for the model ADGP assuming as fiducial 

value D.A = (below) and Ha = 1 — H m (above). 



Model 


constraints 


GDE1 


= 0.022; or(E a ) = 0.008 


GDE2 


cr( 7 ) = 0.036; a(E 2 ) = 0.024; <x(E 3 ) = 0.076 


DGP 


cr( 7 ) = 0.015 


ADGP 


ct(Oy)acdm = 0.051; o-(Ha)dgp = 0.049 



Table II: Constraints on the parameters in terms of the standard deviation a (benchmark survey, fully marginalized). 
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Appendix A: PERTURBATION EQUATIONS 

In this appendix we recapitulate the basics of of first order perturbation theory, both to collect the most important 
equations in one place and to define our notation. We limit ourselves to scalar perturbations, so that we can use the 
Newtonian gauge (also known as longitudinal gauge), following (23|. The metric perturbations are then defined by 
two scalar potentials tp and 0, 




ds 2 = a 2 [- (1 + 2f) dr 2 + (1 - 2<P) dx t dx'] 



(Al) 
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Figure 11: The non linear power spectrum for ACDM and DGP model for one single bin (bin 5). The red (central) solid line 
and the magenta dashed line show the convergence power spectrum for the ACDM and the DGP model respectively. The 
shaded area (delimited by solid blue and brown lines) shows the noise errors on the ACDM convergence power spectrum. 



The total energy momentum-tensor for a perfect fluid is given by: 

= (p + p) u»u v + p g» v (A2) 

where p and p are the density and the pressure of the fluid respectively and the is the four- velocity. We parametrise 
the averaged pressure via an equation of state parameter w, 

p = wp. (A3) 

The perturbed energy-momentum tensor in Newtonian gauge can be written as 

T ° = -(P + Sp) (A4) 
ikT£ = -ikT° = (l+w)p6 (A5) 
Tj = {p + S P )S l J +U i j (A6) 

where 9 can be thought of as the divergence of a velocity field, and IK = T- — S' l jT^/3, the traceless component of the 
space-space part of the energy momentum tensor, represents an anisotropic stress. The scalar part of the anisotropic 
stress a is related to Ilj through 

1 

— ( 

3 



{p + p)a=-(k i k j --6 i j nj (A7) 



We additionally introduce a different velocity variable V = (1 + w)6 which is better behaved at w = — 1 [24j]. In these 
variables the perturbation equations for the fluids are 

Si = 3(1 + --^--(—~ wa) (A8) 



Ha 2 a \ pi 



17 



V- k 2 5n k 2 h 2 

v > = -« " + 5?* + C + «>ffi?* - H*» < A9 > 

where the prime denotes the derivative with respect to the scale factor a. The fluid perturbations are all linked by 
their coupling to the gravitational potentials, 



Pidi + o-—^ 



(A10) 



fc 2 (</>- V) = 12TrGa 2 Y / (p l + Pl )a z (All) 

i 

where the first equation is a combination of the — and — i Einstein equations of [23| . It is often convenient to 
replace the density contrast S by the comoving density perturbation, 

A = 5 + 3— (A12) 

a k z 

and the pressure perturbation Sp is often parametrised with the rest-frame sound speed c 2 , 

5p = c 2 Sp + 3aH(c 2 s -c 2 a )p^ (A13) 

where c 2 = p/p is the adiabatic sound speed. 

As an example, collisionless cold dark matter has zero pressure (w m = 0), vanishing sound speed c 2 m = c 2 a m = 
and no anisotropic stress a m = 0. The perturbation equation for the matter fluid then become: 

V k 2 

V = — + — jV- (A15) 



Appendix B: THE LENSING FISHER MATRIX 



Here we discuss how the lensing Fisher matrix is modified in the general case. Let us briefly recall the main 
equations for weak lensing studies. The convergence weak lensing power spectrum can be written as [il| 



r°° a z 

Pij (£) = Hi ^ — - Wi {z)W 3 (z)P nl 



r{ z y 



(Bl) 



where P n i[Pi(k, z)] is the non-linear matter power spectrum at redshift z obtained correcting the linear power spectrum 
Pi(k, z). In flat space we have: 



-fi m ^(z)(l + z) 



Fi(z) 



dz. 



n t (z s )r(z, z s ) 
r(0,z s ) 

m(z) = Di(z)/ / Di{z')dz' 



H{z) = H E(z) 
dz' 



r(z,z s ) 



E(z') 



(B2) 

(B3) 

(B4) 
(B5) 
(B6) 



where Di(z) is the radial distribution function of galaxies in the z-th z-bin. We assume an overall radial distribution 

^(z) = z 2 exp[-(z/z ) 1 - 5 ] (B7) 

The distributions Di are obtained by binning the overall distribution and convolving with the photo-z distribution 
function. 
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The Fisher matrix for weak lensing is given by 

F aP = f skv ]T (2 ^ +1)A ^ a(^), a C^ 1 fl(J^ T> ),ffC-^ (B8) 

where the cosmological parameters are p a and partial derivatives represent d/dp a , and 

Q fe =P jfe + ^ fe ( 7 L)^ 1 (B9) 

where ji nt is the rms intrinsic shear (we assume (jfnt) 1 ^ 2 — 0-22 [45|]) and 

1 80 

Uj = 3600d( ) 2 h 3 (BIO) 

7T 

is the number of galaxies per steradians belonging to the z-th bin, d being the number of galaxies per square arcminute 
and hi the fraction of sources belonging to the i-th bin. 

As we have seen in the main text, we can parametrize a large number of modified gravity models by the linear 
growth factor and by the combined effect of the modified Poisson equation and the anisotropic stress (the function 
£). So we have that the convergence spectrum can be written as 



f°° dz 

Pij{£) = Hi J -^W l {z)W ] {z)P n 



Q 2 (i + f)W^ 

2 \ r {z) 



(Bll) 



where the functions rj and Q will in general depend on k (and therefore on £) and z. Notice moreover that the matter 
power spectrum P n i depends on the linear growth function which itself is a function of the background expansion 
H{z) and of the functions Q(k,z) and rj{k,z). 
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